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ABSTRACT 



Quantum scattering at zero energy is studied with stochastic methods. A path integral 
representation for the scattering cross section is developed. It is demonstrated that Monte 
Carlo simulation can be used to compare effective potentials which are frequently used in 
! multiple scattering with the exact result. 



1. Introduction 



Multiple scattering off nuclei is in general a complicated many body problem, as target 
and projectile degrees of freedom are strongly coupled. The standard method for treating 
multiple scattering problems is the construction of an effective one-body optical model 
C\l ' potential by eliminating the target degrees of freedom. Optical potential calculations have 

been widely and very succesfully used in the past 0. 
Q\ . Despite their phenomenological success, there are severe shortcomings of these models. 

One example is the spectrum of kaonic atoms, where the shifts of the lowest level require 
a repulsive real part of the optical potential, in contrast to the results of conventional fits 
0. Another problem, which is conceptually even more severe, is the absence of reliable 
methods for calculating inclusive cross sections, for which optical potential models can not 
be applied at all, as they are based on restriction of the target Hilbert space. 

Because of these problems, alternative methods have to be studied. For calculating 
ground state properties of a many body system beyond perturbative or mean field approx- 
imations, stochastic methods are well established. Starting from a path integral expression 
for the density matrix, an algorithm of Metropolis type or Langevin simulation is used for 
path sampling ||. The advantage of these approaches is that they provide results, which 
are in principle exact and can be used to develop better analytical understanding of the 
physical system under investigation. 

2. Path integrals and scattering observables 

We start from a Hamiltonian 

H = Hi nt + J— + V 
2m 
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which can be decomposed into an internal target Hamiltonian i?j„ ( , a projectile kinetic 
energy and a projectile-target interaction V. In case the projectile has no bound state, 
the ground state of the system is the zero projectile momentum scattering wave function 
*&o,k=o{ x 7 <?)• m t ne low temperature limit, the density matrix of the system is dominated 
by this state @: 



(l \ 3/2 

p(x'q, xq\/3) = (x'q'\e^ H \xq) e~^° l^-j * ,k= (x', q'H, k=0 (x, q) (1) 

q denotes the target and x the projectile degrees ao freedom. E is the target ground state 
energy. In pure bound state problems, convergence is controled by the energy gap between 
first excited state and ground state. Here, the ground state of the system lies at the edge 
of a continuum. This fact manifests itself in the /5 3//2 factor in front of eq.([L|). The slow 
convergence, as compared with bound state problems, requires rather long times (3, which 
makes it necessary to choose observables and path sampling techniques carefully. 

From the left hand side of (|l|), a path integral expression can be derived || in the 
standard way. Path sampling methods, however, do not yield the path integral directly, 
but give only paths sampled according to the normalized functional 

P\ x(t ) a(a)] = ^MS[x(t),q(t)}) 

nx[t),q{q)\. j d[x]d[q]eM _ S[x{t):m y W 
This difficulty can be solved by measuring the following functional: 

Q[x(t), q(t)} = exp f dt[V(x(t),q(t)) - U(x(t))] (3) 



The interaction between projectile and target is removed from the numerator and replaced 
by an effective interaction U, which only acts on the projectile degrees of freedom, like 
in conventional treatment of multiple scattering physics. The advantage here is, that 
the stochastic process can be used to test the quality of the effective potential U . The 
expectation value of O in P is in the limit j3 — ■> oo: 

(O) P =°° ® ^°^ k =°( x '^k=o( x ) 



^o,k=o(x'q')^l k=0 {xq) 

where $o is the target ground state and ipk=o is th e projectile scattering wave function to 
the potential U. 

3. Example: Potential Scattering 

To demonstrate the feasibility of this type of calculations I discuss potential scattering. 
In this case numerical integration of the Schroedinger equation provides exact results. 
Scattering off a Gaussian potential V g is considered. As reference potential a square well 
V w is used: 



V g (r) = V exp(-ir 2 ) V w (r) = 3 V™ V ° ' r < ^ 

, r > v5 
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Figure 1: Ratio of cross sections of square well potential and Gaussian potential sr as 
function of potential strength Vq. The line is the exact result, the data points are a 
stochastic calculation at (3 = 100. The path was subdivided into 200 time intervals 
A/3 = 0.5. The projectile mass was set to m = 1. 

V w plays the role of the effective potential U in (^|). The parameters of the square well are 
fitted to reproduce the first two nonvanishing moments of the Gaussian potential. Fig.|] 
shows the ratio of the cross sections of the two potentials sr := oy w jay g as a function of 
Vo. The line is the exact result, the data points are obtained from a stochastic calculation 
at /3 — 100. Path sampling was performed with a simple forward Euler scheme Langevin 
algorithm M. 2.5 x 10 5 paths were used to measure O after an equilibration run of 2.5 x 10 4 
updates. 

In the range of Vq where the reference potential is already a good guess, the stochastic 
calculation reproduces the exact result within 2%. Where this is no longer true, results 
become worse. For Vq = —0.5 the stochastic result deviates about 17%. Note that at this 
value of Vq the cross section is already a = 274, because V g has the first bound state at 
Vq = —0.66. There are two ways to improve the results in this region: One possibility is 
to increase (3. This would require much longer calculation times, as the autocorrelation 
time of the Langevin algorithm increases like f3 2 . The other possibility is to improve the 
reference potential. By adjusting the parameters of the potential successively, one can 
obtain results for the stochastic calculation which do not deviate more than a few percent 
from the exact result, although the same simulation parameters as for the first calculation 
were used. 

4. Discussion 

A new method for calculating elastic cross sections at zero projectile momentum was 
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presented. The crucial point is that this method relies on the comparison of the full 
problem with a reference problem. This makes it possible to study questions related to the 
construction of effective potentials in nuclear multiple scattering by computer simulations, 
which seems to be the natural way for a nonperturbative treatment of many body problems. 
The method can be extended to nonzero momentum by exploiting information contained 
in the (3 dependence of observables. Work in this direction, as well as multiple scattering 
calculations, are currently in progress. 

A severe shortcoming of this method is that at this stage it is not possible to calcu- 
late inelastic or inclusive cross sections. Developement of stochastic methods for these 
problems seems to be promising, as scattering observables will not depend strongly on 
individual nuclear states due to summation over final states. The simple structure of ex- 
perimental data like e.g. energy loss spectra strongly supports this conjecture. 
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